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K-THEORY AND PSEUDOSPECTRA FOR TOPOLOGICAL 

INSULATORS 

TERRY A. LORING 

Abstract. We derive formulas and algorithms for Kitaev’s invariants in the peri¬ 
odic table for topological insulators and superconductors for finite disordered sys¬ 
tems on lattices with boundaries. We find that FC-theory arises as an obstruction 
to perturbing approximately compatible observables into compatible observables. 

We derive formulas in all symmetry classes up to dimension two, and in one 
symmetry class in dimension three, that can be computed with sparse matrix algo¬ 
rithms. We present algorithms in two symmetry classes in 2D and one in 3D and 
provide illustrative studies regarding how these algorithms can detect the scaling 
properties of phase transitions. 


1. Approximately compatible observables 

Compatible observables are given by a rigid definition. If they act on finite Hilbert 
space, the requirement is a basis of vectors that are completely localized for each Aj, 
so XjV = AjV for some scalars. If we repeatedly prepare precisely the same state, 
we can hope to get no variance in whichever Xj we then measure. This sounds more 
like a math theory than any laboratory. What if we repeatedly prepare approxi¬ 
mately the same state and approximately measure one of the Xj in turn, and each 
series of measurements shows little variation? Might this happen because the Xj are 
approximately compatible ? 

In the end, we expect this to mean the commutators [Xj,Xk] = XjX^ — X^Xj 
are small, preferably in the operator norm. However the initial definition should 
involve something like small variance of states. More simply, we can seek common 
approximate eigenvectors. A fundamental quantity would seem to be, for an n-tuple 
of scalars Ai,..., A^, 

(1.1) min max ||AjV — Ajv|| . 

||v||=l j 

If this quantity is small enough, often enough, we could declare these observables to 
be approximately compatible. 

This seems like a nearly impossible minimization, so we seek a proxy. Commuting 
operators have a nice joint spectrum called the Clifford spectrum, so we apply the 
definition of Clifford spectrum to tuples of matrices with relatively small commutators 
and see what happens. The resulting joint spectrum has nice theoretical properties, 
a beautiful relation with A-theory, but remains difficult to compute numerically. 

1 
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Generalizing pseudospectrum to what we call the Clifford pseudospectrum, we find 
an efficiently computable approximation to Equation |1.1[ 


Definition 1.1. Suppose Xi through are Hermitian matrices. Let Li,..., be 
any Hermitian representation of the relations for C£(ifl{C), meaning E* = Fj, Ej = 1 
and EjEfc = —E^Ej for j ^ k. The Clifford e-pseudospectrum of (Xi,..., Xd) is 


A,(Xi,...,X,)= Ae 




-1 


> e 


-1 


with the convention 0“^ = oo and ||<S'“^|| = oo whenever S is singular. The Clif¬ 
ford spectrum of (Xi,..., X^) is Ao(Xi,..., X^), also denoted A(Xi,..., X^). The 
complement of the Clifford spectrum we call the Clifford resolvent set. 


We will use the notation 

H(Xi,...,X,) = 5]X,-®E,- 

and 


BffX,, ...,Xd) = 5(Xi - ... ,X„ - XJ) 

= H(Xi,...,X„)-H(A/,...,A,/)). 


For example, 

A(Xi,..., X,) = {A e M-" I BffXi,. ..,Xd) is singular} . 

The representations of Cidfi{C) are not complicated, so it is routine to show this 
definition does not depend on the choice of the Ej. When we get to X-theory we will 
need to keep the matrix size as small as possible to avoid multiplicity in the spectrum 
of5;,(Xi,...,X,). 


Lemma 1.2. Suppose Xi,..., X^ are Hermitian n-by-n matrices. IfX is in Ae(Xi,..., X^ 
then there is a unit vector v in CC with 


IIVv 


AjVll < 



£“ + EiiiY’W||| 


for all j. 


Proof. Assume we have selected the Ej in Mg(C) where g is minimal, so g 
Without loss of generality, we assume A equals 0, and 

e=||(fi(Xi,...,X,)- 




1-1 



Since B(Xi,..., Xd) is Hermitian, e has the alternate description as the absolute 
value of the smallest eigenvalue of B(Xi,... ,Xd). Let z be a corresponding unit 
eigenvector. Since 


(1.2) H(Xi, ...,Xd)^ = J2 ® ® r,Tfc 
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we make the estimate 


<\\B{X^,...,Xdfz\\ + 


'^[Xj,Xk] 0TjTk ) Z 




Let 


z = 


Zi 


L Zg 


and let r be an index maximizing ||zr||, so we have ||zr|| > 1/g- Let v = z^/ ||zr||- 
Since 




we have 




<9 




Since Xj < X| we find 


(YYv><(E^|v,v)< 




and so 


(A>.v)<j(£" + ^|||V.A-J||V 

V j¥/o / 


□ 


Lemma 1.3. Suppose Xi,... ,Xd are Hermitian n-by-n matrices. If 

\\XjV - Ajv|| < e 

for all j, then A is in Ag/(Xi,..., X^) where 


e = 


d + 1 

\ 


A^ + ElIK'^'dll- 

ii^k 


Proof. Again we use linearity to reduce to the case A = 0. If ||Xjv|| < e for all j then 
let 


z = 


V 
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Using Equation |1.21 we find 








= \ 9 






£ v^EU'MI + ElllY.^i 


3¥^k 


j^k 

This gives a lower bound on the norm of (i?(Xi,..., X^))^, specifically 

V j¥/o / 

Since B{Xi ,..., X^) is Hermitian, we conclude 

||(B(Ai,,,, ‘ < ^ |||A,, AJ||, 

Ji^k 


□ 


Remark 1.4. Lemmas |1.2| and |1.3| tell us that for almost commuting Hermitian ma¬ 
trices Xi,... ,Xd we can get an approximation to the quantity in Equation |1.1| by 
computing 

\\{B^{X^,...,Xa))-^\\~\ 

In a numerical setting, we can compute this easily. Eor example, we can compute the 
absolute value of the eigenvalue of B(Xi ,..., Xa) that is closest to zero. This matrix is 
Hermitian, and typically sparse, so standard algorithms work well for modest matrix 
sizes. The algorithms typically compute an associated (approximate) eigenvalue, so 


we have a way to construct vectors that come close to the minimum in Equation 1.1 


As we push the matrix sizes larger, we will need to do better. Still, estimating the 
norm of an inverse is a fairly standard problem in numerical analysis. One issue is 
that it is hard to differentiate an eigenvalue at zero from one close to zero. This is 
why we turn to the pseudospectrum. If we are computing the function 


||(Ha(Ai,...,X,)) 


-1 


1-1 


we need to set a value e just above zero and regard all values below that as equal. 
This is very reasonable, as we are modeling simultaneous approximate measurement 
when true simultaneous measurement is impossible. 
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Example 1.5. If Xi, ..., Xd are commuting Hermitian matrices then A(Xi, ..., Xd) 
equals the usual joint spectrum. This is an immediate corollary of Lemma |1.2| and 
Lemma 11.31 


Example 1.6. If A and B are Hermitian, then (Ai,A 2 ) is in A{A,B) if and only if 
Ai + iX 2 is in the spectrum of A + iB. So the Clifford spectrum of a pair of Hermitian 
matrices is finite. For positive e we can show that(Ai, A 2 ) is in Ag(H, B) if and only if 
Ai + iX 2 is in the usual pseudospectrum of A + iB. However need the convention 

a,{Y) = {a G C| ||(«-y)"^|| > } 

and not the convention with strict inequality, as in the excellent book 133 by Trefethen 
and Embree. To see the connection, we temporarily use 


ri = 


0 1 
1 0 


, r2 = 


0 i 
-i 0 


so that 




0 (Ai + iA 2 ) — {A + iB) 

((Ai + iX2) — (A + iB)) 0 


Often the better choices here are 


ri = 


0 1 
1 0 


, r2 = 


1 0 
0 -1 


as this keeps 


B) 


B — A2 A — Ai 

A — Ai —B + A 2 


real when A and B are real. Then we are able to produce real joint approximate 
eigenvalues for A and B by finding near null vectors of B(^XiM)i^^ 

Next an example where the Clifford spectrum is an infinite set. For d = 3 the clear 
choice for the Tj is Fi = F 2 = F 3 = so that 


B{X,Y,Z) 


Z X-iY 
X + iY -Z 


as was done in previous work with Hastings ini 

Example 1.7. A nice example, computed by Kisil EH. shows us that the Clifford 
spectrum for three Hermitian matrices is radically different from the Clifford spectrum 
of two Hermitian matrices (as defined below), as it need not be a finite set. We 
compute A(cra;, ay, a^) with help from a symbolic algebra package. The “characteristic 
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polynomial” here is 

det {B{ax — rl, ay — si, Gz — tl)) 



1 — t 

0 

—r -\- is 

0 

\ 


0 

-l-t 

2 

—r -|- is 



—r — is 

2 

— 1 + t 

0 



0 

—r — is 

0 

1 +1 

/ 


= (r^ + + t^ + if - 4. 

This means K{ax,<Jy,<Jz) is the unit sphere. 

While investigating D-branes, Berenstein Malinowski [3] took the preceding exam¬ 
ple further. In that setting, the position observables do not commute. They looked 
at higher spin representations and computed the Clifford spectrum, again a sphere. 
In fact they were interested in a subset of the Clifford spectrum that needs some form 
of i^-theory for its definition. 

Where an index, and eventually i^-theory, arise is easily seen in Example |1.7[ Let 
us examine what is going on at two points in the Clifford resolvent set, the origin and 
(0, 0, 2 ). We find that 

1 0 0 0 ' 

0-120 
0 2-10 
0 0 0 1 

which has a single eigenvalue at —3 and a triple eigenvalue at 1. On the other hand 

-3 0 0 0 ' 

0-320 
0 2 10 

0 0 0 3 

has spectrum 

|-l-\/ 8 ,-l,-l + \/8,3} 

and so the same number of positive and negative eigenvalues. As we vary A the 
eigenvalues move continuously. It follows that for any A inside the unit sphere 
B\{ax,(Jy,<Jz) will have just one positive eigenvalue. For A outside the unit sphere 
B\{ax,(Jy,<Jz) will have exactly two positive eigenvalues. The Clifford resolvent set 
contains information and we will see that from a computation of B\[Xi,... ,Xd) at 
a single point we can make predictions on the size of the Clifford spectrum. This is 
the mathematical essence of bulk-edge correspondence. 

We now require that our Ej are selected in matrices of minimal size. In fact, let us 
consider d = 3 and use the Pauli spin matrices, as above. Recall that for an invertible 
Hermitian matrix Q the signature of Q is the number of positive eigenvalues (with 


-®(0,0,2) W) 


Boiax, Gy, Gz) 
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multiplicity) minus the number of negative eigenvalues of Q, denoted Sig((5). Note 
the signature is always even for even size matrices. 

Definition 1.8. If A is not in A(X, Y,Z) then the index of this triple at A is 

lnd^iX,Y,Z) = ^SigiB^iX,Y,Z)) 
which is in the abelian group Z. 

Remark 1.9. We us primarily the notation of pure mathematics. In particular * 
indicates the conjugate transpose of a matrix, and is a special instance of the * 
operation in a C^-algebra. 

Example 1.10. Consider a finite model of a two-dimensional Chern insulator on 
square lattice. That is, with zero for boundary conditions. The Hamiltonian we 
consider a tight binding model, where there are two orbital types at each site on a 
square lattice. We have creation operators Cm,n,p and Cm,n,s at site (m, n) in band P 
or S. Let Cm,n be the sum of these two types of create at the same site. The periodic 
Hamiltonian is 

^ ^ ^m,n ( W ^y) Cm+l.n T h.C. 

m,n 

+ ^ ^ {~^z ~ '^^y) + h.C. 

m,n 

^ ^ ^m,n ((3 + f^m,n) ^z) 

m,n 

where is drawn with uniform distribution from [—y, y] where N sets the disor¬ 
der level. This is the model used for a Chern insulator as part of the numerical study 
done with Hastings ra. which was essentially the spin-up only part of the model for 
an HgTe quantum wells given in m- If we use lattice position (roughly Angstroms) 
in defining our position operators, we find ||[iL, X]|| and ||[i/, F]|| rather large, about 
6 . We work with the triple {riX,riY,H), although we plot our results using lattice 
units. For this example, r] = 0.5 was selected as a value for which the computed 
approximate eigenvectors we spread out roughly one nanometer in position. We cal¬ 
culated the e-pseudospectrum, and also the index at many positions at the Fermi 
level. For better viewing, the X coordinate of the A was truncated to [—2, 2] and 
energy coordinate to [—2.5, 2.5]. The full energy spectrum is roughly [—7,7]. This 
portion of the pseudospectrum and labeled resolvent is shown in Figures [Tljll.SI with 
an increasingly large random disorder. The pseudospectrum is calculated at 5 grid 
points per unit and e = 0.05. 

2. Almost commuting matrices, the ten-fold way 

We have many potential sources of almost commuting matrices, but now focus on 
situations most relevant to topological insulators. A lattice model of a D-dimensional 
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Figure 1 .2. A Chern insulator on a 18-by-18 lattice with disorder at 4. 

topological or ordinary insulator needs D+1 Hermitian matrices to be described. The 
D position operators Xj will commute with each other and almost commute with 
the Hamiltonian H. In most situations our formulas will work just as well if the Xj 
almost commute, so we do not always require the Xj to exactly commute. 

We now consider the ten symmetry classes in the Atland-Zirnbauer [T] classifica¬ 
tion. Depending on the symmetry class, we may have up to three symmetries. Time 
reversal will be denoted T and be antiunitary and commute with all D + 1 matrices. 
Particle-hole conjugation will also be anti-unitary, will commute with position and 
anticommute with the Hamiltonian. The symmetry S will be unitary, commute with 
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position and anticommute with H. If all three are present, then S is the product of 
the other two, which commute. All symmetries are of order two. 

Following Kitaev | 22 | , we will focus on situations where there is a localized spectral 
gap. A minimal interpretation of this is 


However, we are dealing with approximate measurement, so a better definition is 
generally 


\\B{Xu...,Xn,H)-^\\<S, 

and||[X,-,Xfc]|| <(52and||[A,-,i/]|| < 82 ■ For each choice of and 82 we are specifying 
a potentially useful collection of systems. 

In lower dimensions we can hope to classify such systems up to homotopy and 
identify computable invariants to detect these homotopy classes. As we move to 
higher dimensions we need to allow for stabilization by adding on trivial systems. 

In addition to homotopy questions, we can ask of a given system is close to another 
system that is in the atomic limit, where the Hamiltonian commutes with all the 
position operators. This is then a special case of a mathematical question. Given 
D + 1 almost commuting Hermitian matrices Xi,, X^^i in specific AZ class that 
almost commute, are there nearby commuting Hermitian matrices in the same class 
that commute. To be a serious question this must be posed in a way that is uniform 
for all matrix sizes. That is, in Theorem |2.1| the selected 5 must work for all choices 
for the matrix size n. 

One instance of this, for dimension 2 in class D, asks the following. Given two 
real symmetric and one imaginary antisymmetric matrices that almost commute, 
can these can be uniformly approximated by commuting matrices, again with two 
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Cartan 

class 

r 

c 

s 

Complex 

A 

— 

— 

— 

AIII 

— 

— 

/ 

Real 

AI 

ToT = I 

— 

— 

BDI 

ToT = I 

CoC = I 

/ 

D 

— 

II 

o 

— 

Dill 

1 

II 

o 

II 

o 

/ 

All 

1 

0 

— 

— 

CII 

1 

0 

1 

o 

/ 

C 

— 

o 

II 

1 

— 

Cl 

II 

o 

o 

II 

1 

/ 


Table 1. The symmetry classes. The mark / means this symmetry 
exists, and SHS = —H and SXjS = Xj for a unitary S, with H 
the Hamiltonian and Xj any of the position observables. The mark — 
means this symmetry does not exists. If the other two symmetries exist 
then S = ToC = CoT- All other formulas indicate the symmetry exists, 
with the stated conditions holding, with SHS = —H and SXjS = Xj, 
or HoH = HoH and 1 ~ o Xj = Xj o X, ov C o H = —H o C and 

CoXj = Xj O c. 


Cartan 



Dimension 



class 

0 

1 

2 

3 

4 

5 

6 

7 

Complex 

A 

z 

0 

Z 

0 

z 

0 

z 

0 

AIII 

0 

z 

0 

z 

0 

z 

0 

z 

Real 

AI 

z 

0 

0 

0 

z 

0 

Z2 

Z2 

BDI 

Z2 

z 

0 

0 

0 

z 

0 

Z2 

D 

Z2 

Z2 

z 

0 

0 

0 

z 

0 

Dill 

0 

Z2 

Z2 

z 

0 

0 

0 

z 

All 

z 

0 

Z2 

Z2 

z 

0 

0 

0 

CII 

0 

z 

0 

Z2 

Z2 

z 

0 

0 

C 

0 

0 

z 

0 

Z2 

Z2 

z 

0 

Cl 

0 

0 

0 

z 

0 

Z2 

Z2 

z 


Table 2. The range of the invariants [231 ES|- 
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being real symmetric and being imaginary antisymmetric. The answer is no, with an 
obstruction in Z, as indicated in table [ 2 | 

We tend to prefer describing the symmetries in terms of operations on matrices m 
So we work with the dual operation JJ that is derived from fermionic time-reversal by 

= 7-1 o gw r 


were 


T 


V 


—w 

w 


V 


■ A 

B ■ 

tt 

■ 

-ht ■ 

c 

D 


_^T 

AT 


In block form, 


The first two columns in Table [ 2 ] are unique. These invariants are just a refiec- 
tion of the homotopy classes of such locally-gapped systems. In these columns the 
invariants do not represent obstructions to a system being close to another system 
in the atomic limit. The mathematics behind this statement is nontrivial. It says 
that approximately measuring two incompatible observables simultaneously is very 
different from doing so with three or more. 


Theorem 2.1. (Lin’s Theorem, 1995 ) For any e > 0 , there exists 5 > 0 such that 
whenever n G N and two self-adjoint matrices H and X in M^(C) satisfy ||iL|| < 1 
and ||X|| < 1 and 

\\HX -XH\\ < S, 

there exists a pair of self-adjoint matrices Hi and Xi in M„(C) such that 

\\Hi-H\\<e, ||Xi-X|| <e 

and HiXi = XiHi. 


Lin’s original proof is difficult, so perhaps a better starting point in the litera¬ 
ture is na. However, Ogata [33 adapted Lin’s original proof to work with more than 
two almost commuting Hermitian matrices in a special case involving macroscopic 
observables. These observables are multiparticle averages that avoid the iL-theory in 
columns two and above in Table [ 2 l 


Conjecture 2.2. Lin’s Theorem remains true of we assume that {X, H) is in any 
Atland-Zirnbauer symmetry class and we require that {Xi,Hi) be in the same sym¬ 
metry class. 

Joint work with Sprensen [301 El] proved that this conjecture is valid in classes AI, 
D, All and C. Of course Lin dealt with class A, leaving the conjecture open on the 
five classes involving two antiunitary symmetries. Even in the complex case, research 
into Lin’s theorem is not over, in particualar looking at quantitative versions, as in 

[Eiiini. 
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The only serious consequence of Lin’s theorem we explore here is the following 
remark about ID systems. However, there are expected to be results regarding the 
classification of 2D systems, along the lines of the results in [SI]- 

Remark 2 . 3 . Consider two almost commuting Hermitian operators H and X. The 
Clifford spectrum of this pair will be a perturmation of the Clifford spectrum of 
a commuting pair, which is a finite set. The Clifford pseudo spectrum computed 
for ID systems has typically been a collection of small disconnected regions. The 
computations done have been limited, but it is expected that the pseudo spectrum of 
a ID system look very different from the mutated spheres we see for 2 D systems. 

3. Dimension zero, index formulas 

Many approaches to the iL-theory of topological insulators rely on spectrally bat¬ 
tened Hamiltonian, or equivalently the Fermi projector. For example, see laiiiiisa 
[HI]- One approach that avoids this is the scattering matrix approach muHi. There 
are many numerical issues related to working with matrices with high degeneracy in 
the spectrum, so we prefer to work directly with the full class of invertible Hermitian 
matrices. 

The homotopy classification in dimension zero is rather standard. However the 
formulas for the invariants are not all standard. Our invariants are only designed 
to work for finite models, but it is anticipated that there will be connections with 
indexes defined to work in infinite volume, such as [iniissEii. 

Theorem 3.1. In each of the ten Atland-Zirnbauer symmetry classes, two invertible 
n-by-n Hermitian matrices are homotopie via invertibles in that class if and only if 
the values of an index (in the group 0 or Z or 1 ^ 2 ) for eaeh matrix are equal. This 
index ean be computed in OipH) time. The index for eaeh class is listed below. 

An important consideration is how these invariants can be computed more quickly 
than 0{rT‘) when the matrices are sparse. We will give brief remarks on sparse 
algorithms below. 

3 . 1 . Class A in OD. There are no symmetries, except H = H* in GL„(C). The 
index is 

ISig(ff) 

and it is just a variation on the spectral theorem that this classifies such matrices up 
to homotopy. 

Remark 3 . 2 . The signature can be computed using the LDLT decomposition, which 
finds a unit lower triangular matrix L and a block diagonal matrix D with 1 -by-l and 
2 -by-2 blocks so that H = LDL*. By Sylvester’s law of inertia, Sig {H) = Sig {D) and 
the signature of D can be found in linear time. Since LDLT is 0{rT‘) we are done 
in this case. If H is sparse, there is a readily available sparse version of the LDLT 
algorithm [Zj. 
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3.2. Class AI in OD. We now have H real. We can view that as the added symmetry 
= H. The index is again 

lsig(/f) 

and the algorithms mentioned in §3.1| are available also in the real case. The essential 
fact in proving that this invariant classihes is that H can be factored as H = U*DU 
with D diagonal with decreasing diagonal terms and U being real orthogonal with 
determinant one. 

3.3. Class BDI in OD. The symmetries here can be taken to be = H and 
HT = -TH for H = H* in GL^(C) with 


r = 


I 

0 


0 


That is, 

0 C 
0 

for C an invertible real matrix. The index in this case is 


H = 


sign (det (C)) 

in Z 2 = {±1}. Recall that two real invertible matrices can be path connected if and 
only if their determinants are of the same sign. 

Remark 3.3. There is a noncanonical choice to be made here, specifically a real unitary 
from C” 0 0 i-A 0 0 C. 

Remark 3.4. The sign of the determinant can be computed using the LU decompo¬ 
sition, which finds a lower triangular matrix L and an upper triangular matrix R so 
that C = LR. Then 

sign (det (C)) = n sign(Ljj) Y[ sign{Rjj) 

which avoids the underflow and overflow associated with computing determinants of 
large matrices. Since LU is 0{n^) we are done in this case. If H is sparse, there is 
are readily available sparse versions of the LU algorithm [H]. 

3.4. Class D in OD. The symmetry here can be taken to be = —H for H = H* 
in GL 2 n(C). The eigenvalues of the real, normal matrix iH will come in conjugate 
pairs so its determinant is positive. This makes the Pfaffian real, and our invariant is 

sign (Pf (iH)). 

The homotopy classification can be understood here in terms of the factorization [T8l 
Theorem 9.4] of H as H = UDU* where U is real orthogonal and and D is diagonal. 

Remark 3.5. The sign of the Pfaffian of K = iH can be computed using a decompo¬ 
sition K = LDLU- where D is tridiagonal. If H is banded, which will happen when 
working with derived Hamiltonians based on ID systems, one can use software for 
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PfafRans by Wimmer [3H|- There is an algorithm, but no available software, for the 
more general case of sparse, real skew-symmetric matrices (H]. 

3.5. Class All, OD. We now have H self-dual, so H* = = H. The index is 

jSig(i/) 

due to Kramer’s doubling. The homotopy classihcation is best understood here in 
terms of the factorization of m Theorem 4.6], H = UDU* where [/ is a symplectic 
unitary and D is diagonal with 


4. Dimension one, index formulas 

Although Lin’s theorem is deep, we can avoid it when we have two incompatible 
observables if we are willing to study systems up to homotopy. Given H and X 
Hermitian matrices, one form of our local gap condition 


0 

X + iH 


{X + iH)* 

0 


is invertible 


translates to the condition W = X+iH is invertible. It is easy to deform an invertible 
to a unitary by the path 

Wt = W { W*Wf 

for t in [0,1]. Extracting the Hermitian and anti-Hermitian parts Xt = -\- \ Wt 

and Ht = — |1W we get a path to a system the atomic limit. At all points 

on the path the commutator [Xf , Ht] will remain small if the initial commutator is 
assumed sufficiently small. One can check that this construction respects the various 
symmetries. For example, in class C we have H'^ = —H and X'^ = —X. This implies 
IT** = W*. Functional calculus commutes with ]] so 

(ITt)" = (W*wy IT* = (ITt)". 


Finally 

x! = Iw; + IWt = Xt 

and 

h! = IWt + Iw; = -Hf 

By this homotopy argument, we can simply check that a formula is invariant under 
homotopy and that correctly classihes systems in the atomic limit, at least up to 
homotopy. 
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4.1. Class AIII in ID . We can assume we have H = H* and X = X* in M„(C), 
with X + iH assumed to be invertible. With 


r = 


I 0 
0 -I 


defining a grading, we have HT = —TH and XF = FX. This means that (X + iH) F 
is Hermitian. It is invertible because F is unitary. The index we use here is 

iSig((X + ^/^)F). 

Consider the case where IT = X + iH is unitary, in this class. This means = U 
where a denotes conjugation by F. This symmetry will hold also for /(U) so long as 
/(A) = /(A). We can select a vertical line in the complex plane, near the imaginary 
axis, that misses the spectrum of U and get a homotopy ft from the identity function 
to the function that maps all to the right of the line to 1 and all to the left to —1. 
Thus we have a homotopy to a unitary that is symmetric. Back in the X and H 
picture, we can assume H = 0. Along with the fact that X is even, we are in the 
situation 

■ A 0 

0 B 


X = 


where A and B are Hermitian. Since H = 0 the index is 

lSig(.4)-lSig(B). 

The reason this is the correct index to classify, where it seems we need two indices, 
is that we can bring back nonzero H and use paths such as 


Ut = 


cos(t) 


sin(t) 


sin(t) 


■ cos(t) 


to increase the signature of B by to while decreasing the signature of A by the same 
amount. 


4.2. Class BDI in ID. We can assume H is odd and real and X is even and real, 
where even and odd is determined by the grading operator 


F = 


I 0 
0 -I 


We can use the index from §4.1[ but that turns out to be slower to compute that is 
necessary. 
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Example 4.1. Suppose a is real, between 
example has 


X = 


a 

0 


— 1 and 1, and set /3 = Vl — 

0 ■ 

a 


Our 


and 


H 


0 13 

P 0 


These have the correct symmetries, and 


(X + iH) r 


a —ip 
ip a 


This has signature 0. It is unitarily equivalent to a real matrix, which is not an 
accident. 


We use the unitary Q = ul -\- cuT where ^ \ ~ \ - We can conjugate by a unitary 

matrix without altering the signature, and we compute 


g (X + iH) TQ* = {ujI + cJT) (X + iH) T {UJl + wT) 

- xr + i^ 


and discover a nice index, 

lsig(,Yr + /f) 


which involves now the signature of a real symmetric matrix. See Remark |3.2| on 
computing signature, especially for sparse matrices. 

The proof here that this invariant suffices is almost identical to the argument in 
Section 14.11 


4.3. Class D in ID. Since = I we can assume H^ = —H and X^ = X. Since 
these are Hermitian, we see X is real and H is pure-imaginary, so X -|- iH is real. We 
are assuming it to be invertible. The index here is 

sign (det (X -|- iH)). 

We know the sign of the determinant classifies real orthogonal matrices up to homo- 
topy. We apply this to X -|- iH and extract the needed Hermitian and anti-Hermitian 
parts. 


4.4. Class Dill, ID. We can assume H is imaginary and self-dual, and X is real 
and self-dual, but in this case we find that even and odd are to be determined by the 
grading operator 


0 -il 
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This ensures that the transpose is conjugated to the dual. We use the unitary Q = 
ujI + cJT where ^ ~ 2 ’ compute 

g (X + iH) rg=^ = + ^r) (x + iH) r {ui + wT*) 

= {ljI + hJT) (X + iH) {uJT + ljI) 

^-'-(x + iH) + ir(x + iH) + i(x + iH)r+'-r(x-iH)r 

= XT+ H 


and we check its symmetries 

(XT + H)* = (r)X + H = (XT + H) 

and since T and H are imaginary and X real, this is imaginary, and so skew-symmetric. 
We can compute a Pfaffian, as before, and our index is 

sign(Pf {i {XT+ H))). 


Here is a sketch of why this invariant classihes. As in class AIII we are able to use 
functional calculus to reduce first to the case of X -|- iH being unitary and then also 
to where H = 0. If X is real symmetric and self-dual and unitary, it can be shown 
that it will factor as 


Q* 


D 


D 


Q 


with g real orthogonal and symplectic and 


1 



-1 


To see how to finish the classification, notice the invariant comes out differently on 


X^ = 


1 


since 


On the other hand 


pf (zx±r) = Pf 


x^ = 


-1 


0 ±1 
Tl 0 


-1 


= ± 1 . 


cos(t) 0 0 

0 cos(t) — sin(t) 
0 — sin(t) cos(t) 

sin(t) 0 0 


sin(t) 

0 

0 

cos(t) 


is a real, self-dual, Hermitian unitary path from I 4 to —I 4 . 
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4.5. Class CII, ID. This is much like the BDI situation. On M 4 „(C) the operation 
corresponding to particle-charge conjugation is 


The grading operator is 


' A By 

' A# 

C# 1 

CD 





The operations corresponding to time-reversal is 

A B y _ \ I 0 A By \ I 0 

C D \ “[o -I \ [C D \ [o -I 


I 0 1 r (7“ 

0 -/ J [ 5# 


/ 

0 



A» -C# ■ 
- £>» 


To check these are the correct type, we notice 

'A By _ _ \ z oita b y \ Z o' 

C D \ ~ [ 0 Z \[c D \ [ 0 Z 


and 

' A B y _ _ \ Z 0 ] \ A By \ Z o' 

C D \ ~ [ 0 -Z\[C D \ [ 0 -Z \ ' 

and since 

z 0 y _ \ z ol^__r/o 

0 Z \ ~ [ 0 -Z \ ~ [O I 

we have C oC = —I and 7~o 7~ = —I- 

Our matrices are Hermitian X and H with X even and X'^ = X, and with H odd 
and H'^ = H. The index we use is 

lsig(,Yr + /f). 

The reason for the extra factor of one-half will evident from the symmetries here. We 
note 

{xr + H)* = rx + H = xr + H 

and 

{XT + Hf = T^X + H = TX + H = XT + H 

so this matrix is “self-dual” and Hermitian, so has Kramer’s doubling. 

Arguing as before, we reduce to the case H = 0 and X unitary. This means 
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with Y and Z both self-dual, Hermitian and unitary. Since 


jSig (AT) = tsig u' isig (A) - tsig (Z) 


we initially seem to be short by one invariant. The path 

cos(t) sin(t) 


cos(t) 

sin(t) — cos(t) 

sin(t) 


sin(t) 


cos(t) 


illustrates how to use nonzero Ht to increase one signature by four while decreasing 
the other by four. 


5. Dimension two, index formulas 


What we are after here are invariants that, for two dimensional finite systems on a 
square, can be quickly computed and that can explain the robustness of gapless edge 
modes in the face of disorder. We only consider disorder that respects the symmetry 
class. However, the simplicity of these invariants should mean they function well for 
disorder that is only approximately invariant under the needed symmetries. 

5.1. Class A in 2D. We have no reality condition and use the index from Section[^ 
at the origin, so 



Notice that if consider 

H' = X <S> (Tx + Y (Ty + H <S> Cz 
as a derived Hamiltonian, it constitutes a class A system in OD. 

5.2. Class D in 2D. We can assume, after perhaps a unitary change of basis, that 
H is imaginary while X and Y are real, all being Hermitian. We select our Tj so that 
H is tensored with ay, which is imaginary. Therefore 


A Y -iH 
Y + iH X 


H' = X Z> az + Y Z> ax + H Z> ay = 


defines a OD system in class AI. That is, it is real symmetric and, by the local gap 
assumption, invertible. In terms of Table [2] this moves us two steps up and two to 
the left. This is similar to the scattering matrix approach of Fulga et al. [H] which, 
in slightly different geometry, moves one step up and one step left. 

We utilize the invariant from Class AI in OD and use 



This is invariant under homotopy, is additive with respect to direct sums, and is 
trivial on trivial systems. It could be the trivial invariant. Perhaps the best way to 
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show these invariants nontrivial is to use them in an numerical study. We do that in 
some cases. Here we derive the existence of a nontrivial example mathematically. 
The standard example pin] of three almost commuting Hermitian matrices 


W. r X Y-iH 
2®^® YYiH X 


±1 


has one matrix, say H, imaginary and the others real. All we are missing is having X 
and Y commute. By the class AI version of Lin’s Theorem we can modify these 
a little to produce Xi and Yi that are commuting real orthogonal matrices. If we use 
a large enough matrix size, B(X,Y, H) will have a large spectral gap, meaning the 
index will be unchanged when applied to (Ai, li, H). 


5.3. Class Dill in 2D. We are given symmetries H'^ = H, X'^ = X and F** = Ton 
top of knowing X and Y are real and H is imaginary, as in Section |5.2[ As we did 
there, we set 

H' = X ® az + Y ® + H (S> (Jy 

This is real, so = H' , while i/d®# = [{' With these symmetry operations, the 

grading operator is Z ® Z. We need to select partial isometries into the 1 and — 1 
eigenspaces for this grading operator, and choose 


■ I 

0 

0 

I 

TX. 1 

I 

0 

0 ■ 
I 

0 

-I 

w = — 

’ V2 

0 

I 

I 

0 


-I 

0 



Now we can use the index from Class BDI in dimension zero and define our index as 


sign (det {WIH'WJ)) . 

We want to see this is not a trivial index. 

Suppose (X,Y,H) is any example from class D, which can exist with both odd 
index and even. Let our class Dill example be the 2n-by-2n matrices 


(5.1) 

Xi = 

■ X 

X 

, n = 

1 1 

1_1 

, H,= 

■ H 

H 



The index is then the sign 

of the determinant of 











X 

0 Y-iH 

0 



I 

0 ■ 

1 

■ / 0 0 

I ■ 



0 

X 

0 YYiH 



0 

I 

2 

0 / -/ 0 


Y + iH 

0 

-X 

0 



0 

I 






0 Y-iH 

0 

-X 



-I 

0 


X Y-iH 
-Y-iH X 



X 

Y-iH' 


Y + iH 

-X 
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SO the 

Remark 5.1. Notice that in this example, the index comes out 1 when we start with 
X and Y and H all commuting. An easy homotopy argument shows we get trivial 
index of 1 whenever we start in class Dill with all three matrices commuting. This 
means we made valid choices for W±. 

5.4. Class All in 2D. We have one symmetry = H, X^ = X and T** = Y. We 
set 

H' = X ® + Y ® ay + H ® 

and find = —H'. We are in a nonstandard version of class D, dimension zero. 

The unitary matrix 

r / 0 0 -il' 

J_ r / -iZ 1 _ J_ 0 I il 0 

^ ~ y/2 'iZ I ~ 0 il 1 0 

[ -i/ 0 0 1 

brings us to the standard picture and our invariant is 

sign(Pf(zg*//'g)). 

We conduct a numerical study that provides ample evidence that this is not a trivial 
invariant. 

Remark 5.2. There are many choices here, including conventions as to the definition 
of the Pfaffian, so that it is easy to program this wrong. Done correctly, the index is 
1 when H and X and Y all commute. 

5.5. Class C in 2D. We have one symmetry = —H, X"^ = X and Y"^ = Y. We 
set 

H' = X <Si az + Y Zi ax + H Zi ay 

so that nonstandard version of class All, dimension zero. 

The matrix 

■ / 0 0 0 ■ 

0 0/0 
0/00 
0 0 0 / 

converts is to the standard dual operation. Our invariant is 

jSig(OJf'(J) = isig(ff')- 

To show nontrivial values of this index are possible, consider any example from class 
D. Let our class C example be 

(5.2) 
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Then 


QH'Q = Q 


Y 


A 

0 

Y -iH 

0 

0 

A 

0 

Y -iH 

Y + iH 

0 

-A 

0 

0 

Y + iH 0 

-A 

A Y 

-iH 


” 

+ iH 

-X 





A Y 

-iH 



Y + iH 

-X 


Q 


SO 

mdciX,,Y,,H,)=mdn{X,Y,H). 

6. Dimension three, index formulas 

We believe in all five interesting classes in three dimensions we can move up two 
and left two in Table and arrive at a useful index. For now, we study this just 
in Class AIL In this class, we have examples from physics research (with Hastings 
CSl) to show this invariant is interesting, and theorems about the iF-theory of real 
C^-algebras (with Boersema [1]) that allow is to identify the invariant. 

The higher dimensional invariants in classes just one antiunitary symmetry, can 
be explained using the theory of real (ungraded) C'*-algebras. This will be discussed 
elsewhere. 

6.1. Class All in 3D. We have one symmetry H'^ = H,X‘^ — X, Y'^ = Y, and 
Z# = Z. We set 

H' = X (S> (Jx + Y ® Uy -\- Z ® Gz 

and 

X' = H®I 

and find and find and = X' so we have a derived ID system in 

class D. The unitary matrix 

■ / 0 0 -i/ ■ 

J_ 0 I il 0 

^ ^ 0 / 0 

[-il 0 0 I 

brings us to the standard picture and our invariant is 

sign (det {Q* {X' + iH') Q)) = sign (det {X' + iH')). 

Again we offer a numerical study that provides ample evidence that this is not the 
trivial invariant. 
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Remark 6.1. The advantage of the left form of the invariant is that an LU factorization 
will be faster and take less memory. 

7. Bulk-edge correspondence 

We now prove a relation between the index values in the pseudoresolvant and the 
pseudospecturm. Essentially, between two points in energy-position space where the 
index changes, there must be a point in the pseudospectrum. So at that point there 
must be a vector approximately localized in position and energy. For weak disorder, 
this will mean a vector localized at the edge and localized at the Fermi level. For 
strong disoroder, things get messier, with the “edge modes” forming a ring around 
the sample, sometimes moved in from the edge. 

Unlike in [12] and [20], for example, we do not have a separate Hamiltonian for 
the edge states. Rather we are classifying the approximate zero modes of the Hamil¬ 
tonian. In a system with weak disorder and nontrivial invariant in the center, these 
approximate zero modes cannot localize in the bulk, but they can and do localize near 
the edge. Since the value of the invariant is robust against disorder, at least disorder 
with the correct symmetry, whatever is happening at the edge is robust. 

We are not claiming the converse. Most likely, in higher dimensions there are robust 
edge effects that can exist without these particular invariants being nonzero. Indeed, 
see the discussion of stabilization in 77-theory in 

We will use ind(Xi,..., Xjj^i) generically for any of the indices described on Sec¬ 
tions |3ll6l We do not assume in this section that the hrst D matrices commute. Our 
convention is to call Xjj^i the Hamiltonian, even if these matrices are not related to 
quantum systems. 

We will assume the needed symmetries on the Fj to correspond the the choices 
made in the definition of a specihc index. For example when D = 3 and we are in 
Class AH, the index is 



/ 

0 

0 

H + iZ 

iX + Y' 

\ 

\ 

det 

Q* 

0 

0 

iX-Y 

H-iZ 



H-iZ 

-iX-Y 

0 

0 



[ 

-iX + Y 

H + iZ 

0 

0 

) 

/ 


so our choices for the Fi, F 2 , F 3 , F 4 are 


-1 

0 

0 

0 

i 


■ 0 

0 

0 

1 ■ 


0 

0 

i 

1 

0 


■ 0 

0 

1 

1 

0 

0 

0 

i 

0 


0 

0 

-1 

0 


0 

0 

0 

—i 


0 

0 

0 

1 

0 

—i 

0 

0 


0 

-1 

0 

0 

•) 

—i 

0 

0 

0 

•) 

1 

0 

0 

0 

—i 

0 

0 

0 


1 

0 

0 

0 


0 

i 

0 

0 


0 

1 

0 

0 


What is essential is that the index can be computed continuously from 


BiX^,...,XD+l) 

whenever this is invertible. This means that given a continuous path 

{Xi^t, ■ ■ ■, XD+l,t) 
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the only way for to change is for it to be undefined at some point because 
is singular. 

Definition 7.1. Assume (Ai,..., are in some symmetry class. If that sym¬ 

metry class does not have any symmetry that anticommute with the Hamiltonian, 
then for any {D -\- l)-tuple A not in A(Ai,..., define 

indx(Ai,..., Ad+i) = ind (Ai — Ai/,..., — \d+iI) ■ 

In the other classes, we do the same except only define this index at a point A when 
Ad+1 = 0. 

Lemma 7.2. If (Ai,... , A£)+i) and (li,..., Y^^i) are tuples of Hermitian matrices 
of the same size then 

ii(ij(A'i,..., AD+i))-‘ir‘ - ii(B(r,..., u+i))“‘ir‘| < 

If (Ai,..., A/j+i) and {Yi,, Ld+i) are in the same symmetry class and 

ind(Ai,...,Ap+i)^ind(F,...,yD+i) 

then 

||(B(A'i,.... AB+i))-‘|r' + ||(B(r... ,yi,+i))-'ir‘| < 

and somewhere on the line segment (Ai^j,..., A^+i^t) between these pairs is a point 
where 

5(Ai,„...,A^+i,,) 

is singular. 

Proof. We apply Weyl’s estimate on the spectral variation of Hermitian matrices to 

B{X,,t,...,XD+i,t) 

at various places along the line defined by 

A,,, = (1 - t)X, + Y,. 

Let’s let a be the eigenvalue of smallest magnitude for B{Xi,... ,Xd+i), so 

a = ±||(H(Ai,...,Az,+i))-'||“' 
and similarly have an eigenvalue b for B(Y^^ ^+i) with 

6 = ±||(5(W,...,F^+i))-'fL 

Without loss of generality, |a| < |6|. The closest eigenvalue to a in the spectrum of 
B{Yi ,..., Td+i) is af least at a distance of |6| — |a| away. This is then a lower bound 
on the best overall spectral pairing, and that in turn a lower bound on the distance 
between the matrices. 
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If the index varies, then at some t the matrix 

is singular, since the index is continuous where defined. We can apply the first 
statement in the result to B{Xi ,..., Xjjj^i) and B{Xi t ,..., j), and then again to 

B{Yi ,..., Yd+i) and , XD+i,t)- These points are collinear so the estimates 

add. □ 

Lemma 7.3. For any two {D + 1)-tuples A and fj, of scalars, 

\\B (Ai/,..., Ad+i/) — B {fill,... ,/iD+i/)|| = d(A, fi) 

where d denotes Euclidean distance. 


Proof. We first notice 

B (All, ..., Xd+iI) — B (pi/, ..., pd+iI) = B ((Ai — pi) /,..., (A/^+i — /iz^+i) I) 
so we need only compute 

\\B (fill, . . . ■ 


By Equation 


1.2 




SO, by the spectral mapping theorem, 




Therefore 

||B(7,/,...,7d+i/)II = 

□ 


So now we know that points in the pseudospectrum with different index must be 
separated by certain distance. We also have a means to estimate how much disorder 
is needed to change an index. The reason nonzero index is thus important, when it 
is found, is that trivial index prevails away from the origin. 

Lemma 7.4. If {Xi,... ,W£)+i) is in any symmetry class, and if 

\\\>\\B(X,,...,XoW\\ 

then A is in the pseudoresolvent and 

ind^ (Xi, ..., Xd+i) = 0. 

Proof. Notice first that B (AiJ,..., X^+iI) has spectrum {± |A|} with trivial index, 
in any symmetry class. We then consider the path B\{tXi,. .. ,tX£)+i) which has 
length 

I|b(V,.....Yd+,)II- 

For the signature along this path to change, it will need to be at least as long as 
|A|. □ 
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We now prove bulk-edge correspondence. To accommodate highly disordered sys¬ 
tems, we allow for a generous interpretation of edge mode. A nontrivial index in 
the center will need to be surrounded, but not too closely, by a ring of approximate 
modes approximately at the Fermi level (assumed here to be zero). We can, for less 
disordered systems, compute the index well away from the center and see that the 
edge modes are really near the edge. 


Theorem 7.5. Suppose (Ai, ..., X^, H) is in any symmetry class and that 

md{X,,...,XD,H) 

is nontrivial and 

R=\\{B{X,,...,Xd,H))-Y" 

is positive. Then for every unit D-tuple («i,..., at some point along the ray 

t ^ \t = {tai ,..., to/5, 0) 

there is a Xg which is in the pseudospectrum. If Xt is in the pseudospectrum then 


R < |A,| < \\Hf + Y, II Alt + E ll[Y.-H|ll •• 


Proof. Suppose Xt is in the pseudospectrum. The claim |Ai| > R follows directly from 


Lemmas |7.2| and |7.3[ We know from Lemma |7.4| that Xt is in the pseudoresolvent 
whenever 

|A/| > ||S(Ai,...,A/5,77)11 . 

We now use Equation |1.2| to get a bound on this norm, 

\\B{Xu ...,Xd, H)\f = ||S(Ai,..., A/5, HfW 


< 


w+YP + EiA.-ff] 


< llfff+ EllYlt + ElllY'-»lll 

If Xt is in the pseudospectrum then 

|A,t < ||ffIt + E IIAlt + E NIY.-fflll ■ 


Again using Lemma 7_A we see that somewhere on this ray the index is trivial. At 
the start of the ray the index was assumed to be nontrivial, so Lemma |7.2| tells us 
that at least one point on the ray is in the pseudospectrum. □ 


Remark 7.6. In the symmetry classes without particle-hole conjugation, we can con¬ 
sider rays that slant up or down in the energy direction. In those classes, the edge 
modes persist above and below the Fermi level. So when 77 = 2 in a class like All, 
the pseudospectrum contains a sphere that surrounds a hole. So starting with lattice 
geometry of a square, we have produced a modihed sphere and nontrivial topology. 
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8. Local and global invariants 


Suppose we have increasing system sizes, say {Xk,Yk, Hk) with consistent spacial 
units, such as nanometers. We then must select rescaling of units % > 0 for all k and 
so work with the observables (riXk,riYk, Hk). Another view point is that we will be 
basing our pseudospectrum and index on 

B(rikXk,rikY, H) — B{r}kXI, A 3 /). 


In terms of approximate modes, this means we are producing v that more or less 
minimize 

max {rjk \\XkV - Aiv|| , rjk ||Wv - A 2 v|| , ||i7fcV - AavH) 
which, for points in the pseudospectrum, will be on the order of 

yXrkV\\[Hk,Xk]h\\[Hk,Yk]\\. 


he. 


max 


1 


\\XkV - Aiv|| , ^ ||FfcV - Aavll , —— \\HkV - AgvH ) C 

yVk 


If we are trying to model a local probe, we would then want to keep ijk constant. 
Ideally it would be set to correspond to modes localized in energy so most energy 
spectrum can fit in the band gap of the clean periodic Hamiltonian and with spread 
in position of a few nanometers. This will to correspond to the expected spacial 
resolution of a scanning tunneling spectroscope. If we hold rjk constant we will have 
a local index and local pseudospectrum. 


If we instead think bulk-edge correspondence, as in Theorem 7.5, then we want 


Bk \\(^B(^fjkXk, f]kYk, Hk^^ 

Vk 


-1 


1-1 


to grow proportionate to ||A7fc|| which, for simplicity, we take equal to ||yfc||. At the 
very least, we want this to grow, so we need rjk \ 0. In keeping with how we rescale 
the periodic observables in the torus geometry, we use % = 770 / ||A^fc|| and call the 
resulting index a global index. 


9. Dimension two, numerics 

We present examples in two symmetry classes, AI and AIL We look at examples 
illustrating the local invariant, and do a study of the disorder-averaged global invari¬ 
ant, looking for the type of phase transition out of topological insulator caused by 
doping. We use Hamiltonians previously used in such studies so we can test the new 
algorithms. 

We have collected ample evidence that the new indices correspond to old indices, in 
three combinations of dimension and symmetry classes, in the special case of modest 
disorder strength and where the Fermi level is in the middle of the gap of the peri¬ 
odic Hamiltonian of clean system. We also get equality for trivial systems, meaning 
systems in the atomic limit. What can we expect to prove here? 
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Figure 9.1. Chern insulator with zero disorder, looking down at the 
pseudospectrum and labled pseudoresolvant at the Fermi level. Panel 
(a) is the 0.1-pseudospectrum of the Hamiltonian with periodic bound¬ 
ary conditions. Panel (b) is the 0.1-pseudospectrum of the Hamiltonian 
with zero boundary conditions. Panel (c) is the Clifford pseudospec¬ 
trum with K-theory labels in the gaps, but only at the Fermi level of 
zero. 


We are discussing invariants that take discrete values on individual hnite systems. 
These are closely related to the it'-theory obstructions to fixing approximate matrix 
representations of C'*-algebra relations to be exact representations [9] . Based on that 
older research, and the more recent work |28] . we expect to soon hnd a rigorous 
proof that the new index formulas agree with some established index in the case of 
weak disorder and with the Fermi level in the middle of a big bulk gap. For the 
more interesting situation, exploring transitions between topological and ordinary 
insulating states, we expect at best a probablistic result. Moreover different indices 
will mark the transitions differently. We hope to have established a connection with 
these new indices and edge modes, and the practicality of working with systems 
of nontrivial size. It will take time to do numerical studies contrasting the phase 
transitions as found by various index algorithms, including [321131] and [I11[I5|. 


9.1. Class AI in 2D. We look again at the model Chern insulator as in Exam¬ 


ple 1.10 We create the usual Hamiltonian H^er with periodic boundary conditions 


as well as H, the Hamiltonian with zero at the square boundary. We look now at the 
e-pseudospectrum only at energy zero and compute the local index, using r} = 0.02 
and e = 0.05. For comparison we compute the spectrum of both H and Hp^r- In 
fact we plot the e-pseudospectrum for small e as this is known to close to the actual 
spectrum and is much faster to compute. 

At the 


We look at an 18-by-18 lattice with increasing disorder, in Figures 9.1 9.4 


high value of disorder substantially, as in Figure [9^ the zero modes have, in places. 

















K-THEORY AND PSEUDOSPECTRA FOR TOPOLOGICAL INSULATORS 


29 




Figure 9 .2. Chern insulator with disorder set at 7. Panel (a) is the 
0.1-pseudospectrum of the Hamiltonian with periodic boundary condi¬ 
tions. Panel (b) is the 0.1-pseudospectrum of the Hamiltonian with zero 
boundary conditions. Panel (c) is the Clifford pseudospectrum with K- 
theory labels in the gaps. The bulk gap is just closing at this disorder 
strength. The small dots (black) are in the pseudospectrum. Some of 
the grid points in the pseudoresolvant are labeled by there index, with 
crosses (red) indicating index 1 and with circles (blue) indicating index 
- 1 . 



Figure 9.3. Chern insulator with disorder set to 8 (left) or 9 (right). 


moved in substantially from the edge. On a 50-by-50, in Figure 9.5, we see better 


how the center of the sample has roughly circular patches of the wrong index, while 
the corner of the sample still has some semblance of a boundary effect. 
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Figure 9.4. Chern insulator with disorder set to 10 (left) or 11 (right). 


10 

5 

0 

-5 

10 

10 

Figure 9.5. Chern insulator with disorder set to 11. Now on an 50- 
by-50 lattice, just in the middle of the bulk (left) and at the top-right 
corner (right). 



Now we look at the global invariant with rj — 4/L when the model is on an L-by-L 
lattice. Holding disorder hxed at 8 we compute the index at the center, but with the 
energy level moving between —5 and 0. The results are shown in Figure |9?7| 

We did this study with dense matrix methods and the formula for the Bott index, in 
joint work with Hastings [22] • We re-ran this study in order to compute the following 
proxy for the inverse of spread of the Fermi projector P = Pef, 


U,P 


f] 


V,P 


(9,1) 


concentration (Pg^) = 1/ 


2 
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Figure 9.6. Chern Insulator by the old Bott index. The number N of 
samples for an L-by-L lattice used in these averages was: L = 20, N = 
10938; L = 30, N = 2634; L = 40, = 1040; L = 60, N = 150. The 

left panel shows the disorder-averaged Bott index with various Fermi 
levels. The right panel shows the disorder-average of the concentration 
of the Fermi projector, as defined in Equation |9.1[ 




Figure 9.7. Chern Insulator by the new index. The number N of 
samples for an L-by-L lattice used in these averages was: L = 20, 
N = 15929; L = 30, N = 21326; L = iO, N = 16853; L = 60, 
N = 10224; L = 80, N = 8659; L = 100, N = 2845. The left panel 
shows the disorder-averaged new index with various Fermi levels. The 
right panel shows the disorder-average of gap localized at the origin, as 
defined in Equation |9.2[ 


where U and V are the unitary operators corresponding to ’’periodic observables” of 
position on the torus. In terms of the usual position observables X and Y, where we 
































K-THEORY AND PSEUDOSPECTRA FOR TOPOLOGICAL INSULATORS 


32 


have 


^ < — 


we can define these commuting unitary matrices as 


1 


U = e- 




V = e' 


lY 


The quantity in Equation |9.1| is expected to be similar to 
(9.2) localGap(F) = \\{B{rjX,rjY,H - 

in the new method. It is slow to compute because computing the Fermi projector is 
slow. 

By the Fermi projector we mean the spectral subspace of i/per corresponding to 
(—oo, Ef\. In the event of a large spectral gap we can prove that P will have relatively 
small commutator with U and V, as the indicator function can be calculated as 
/(i/per) for / with reasonable Fourier transform. There is also what is called a 
mobility gap [30] , where there is a region of the spectrum around the Fermi level that 
is filled with eigenvalues that have well localized eigenstates. In that case as well, the 
commutators [P, U] and [P, V] tend to be small. 

There is an integer we can calculate here, the Bott index. Let U = PUP + [I — P) 
and V = PVP + (I — P) and define the index 




( Trace ( — log {VUV*U*) 
\ yziri 


which can be proven to be an integer. This integer will be zero when U, V and P are 
close to a commuting triple of matrices and when it is nonzero such as approximation 
is precluded. Given a full eigensolve of Hp^r, if assemble all the low-energy eigenstates 
to form a non-square matrix W with WW* = P, then a reformulation of this formula 
is 

SR ^Trace log (w*V PU PV* PU*W 

We can compute this from just the eigenvalues of W*V PU PV* PU*W . We do need 
up to half of the eigenvectors of iLper, so an algorithm for the Bott index is easy to 
implement in O(n^) time, but no better. 

Both methods show a 


The results using the old method are shown in Figure 9.6 


sharpening transition, with the new algorithm able to work with large lattices. 


9 .2. Class All in 2D. Now we look at the model used in the the numerical study 
done with Hastings ra. which was the model for an HgTe quantum well given in 
ra. We keep the same disorder and the same strength of the Pbia term that breaks 
in version symmetry as in the old study ESI- 

We are just claiming proof of concept, the our formula in a possible replacement for 
the Pfaffian-Bott index. We can’t get to much larger matrices than before because 
our algorithm to compute the sign of the Pfaffian uses dense matrices. We hope 
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Figure 9.8. Class All insulator with disorder at 0. Panel (a) is the 
O.l-pseudospectrum of the Hamiltonian with periodic boundary condi¬ 
tions. Panel (b) is the 0.1-pseudo-spectrum of the Hamiltonian with 
zero boundary conditions. Panel (c) is the Clifford pseudospecturm 
with A'-theory labels in the gaps. 


this data will inspire the production new software implementing the sparse matrix 
factorization in [^. 

We look again at the local index and pseudospectrum, with r] = 0.5. Figures |9^ 


9.11 show this with increasing disorder. The red x indicate index —1 while the blue 


o indicate index 1. 

Now the global index. 


Here we set rj = 4/L for an L-by-L lattice. Figure 9.12 


is reproduced from m- We cannot work with larger matrices in this symmetry 
class because we do not have software to compute the needed factorization of sparse 
antisymmetric matrices. The results of the new formulas are shown in Figure |9.13 


We see that the transition is probably not sharpening as system size increases, but it 
is hard to tell without the larger system sizes. Notice that 40 lattice units is roughly 
two nanometers. If we are modeling films of roughly one nanometer thickness, we 
ought to be looking at L ~ 200. Such a size may be in reach of the sparse algorithm 
as soon as that is available. 


10. Dimension three, numerics 

We present examples in one symmetry class, AH. 

10.1. Class All in 3D. The first numerical study ra in 3D of the effect of doping a 
topological insulator used an index that worked for periodic boundary conditions. It 
also involved calculating the sign of a determinant, but involving polynomials in three 
variables that approximate a degree-one mapping of a three-torus to a three-sphere. 
The geometry here forces these polynomials to have degree eleven, and the resulting 
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Figure 9.9. Class All insulator with disorder set to 7. Panel (a) is the 
O.l-pseudospectrum of the Hamiltonian with periodic boundary condi¬ 
tions. Panel (b) is the 0.1-pseudo-spectrum of the Hamiltonian with 
zero boundary conditions. Panel (c) is the Clifford pseudospectrum. 
At this strength of disorder, the bulk gap is just closing. 



IQ 



-IQ I-'-■-'-' 

-IQ -5 Q 5 IQ 


Figure 9.10. Class AH insulator with disorder at 8 (left) and 9 (right). 


algorithm was slow. In contrast, as we are working with a cube not a three-torus, we 
have a substantially faster algorithm. 

The local index was defined using rj = 0.25. The results are shown in Figures 
10.11110.41 The red X indicate index —1 while the blue o indicate index 1. 
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Figure 9.11. Class All insulator with disorder at 10 (left) and 11 
(right). 



Figure 9.12. Phase transition in 2D. The larger plot shows the disor¬ 
der averaged PfalRan-Bott index. Reproduced from 


The global index was defined using 77 = 4/L for an L-by-L-by-L lattice. Figure 7.2 
of [18] is replicated here as Figure 10.5, This data was sufficiently noisy that no real 
conclnsions abont scaling could be made. The left panel in Figure [10^ shows how we 
can generate much cleaner data with the new algorithm by using more samples and 
larger systems. It appears now the transition from this 3D topological insulator to 
an ordinary insulator is not sharp, but larger system sizes are needed to clarify this. 
These can be studied on existing machines, but the processing time needed will be 
significant. 
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Figure 9.13. Phase transition in 2D. The left panel shows the disorder 
averaged new class All index. The right panel shows the disorder- 
average of gap localized at the origin. The number N of samples for an 
L-hj-L lattice used in these averages was: L = 20, N = 3544; L = 30, 
N = 2651; L = 40, = 481. 



Figure 10.1. Class All insulator in 3D with no disorder. Showing a 
slice at z = 0 and at the Fermi level. Using a 9-by-9-by-9 lattice. 


11. The algorithms 

Matlab code, with instructions on how to produce many of the figures in this paper, 
will be made available at a data repositorjQ For the larger system sizes, the study of 
the global index was done on multiple computing nodes, each with 32 cores and 64GB 
of random access memory. However the local index, at the system sizes illustrated in 
the figures, can be explored using less than a day on a desktop with 4 cores and only 
8GB of random access memory. 


^ https://repository.unm.edu/hanclle/I928/23449 
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Figure 10.2. Class All insulator in 3D with disorder at 11. Showing 
a slice at z = 0 and at the Fermi level. Using a 9-by-9-by-9 lattice. The 
bulk gap is just closing at this disorder level. 




Figure 10.3. Class All insulator in 3D with disorder at 12 (left) and 
13 (right). 


A good example here is the algorithm for the global class All invariant in 3D. The 
formula in equation |7.1| tells us we need to first compute a matrix 


A = Q* 


0 0 H + iZ iX + Y 

0 0 iX-Y H-iZ 

H -iZ -iX -TO 0 

-iX + Y H + iZ 0 0 


that is real and sparse. We need to compute its spectral gap 


Q 


-1 
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Figure 10.4. Class All insulator in 3D with disorder at 14 (left) and 
15 (right). 



Figure 10.5. Phase transition in 3D, as computed via the old method. 

This figure is replicated from [T^. Shown is a plot of average index for 
L = 6, 8, 10, 12 with L x L x L lattices. Each data point is an average 
of 1700, 1400, 600, 400 samples, respectively. 

as well as the sign of its determinant. 

We rely on LU algorithm |6], as implemented in Matlab, to factor A as 

(11.1) A = R-^P*LUQ* 

where R is diagonal, P and Q are permutation matrices, L is lower triangnlar, sparse 
with nnit diagonal, and U is upper triangular and sparse. The determinant of L is 
one, and so 

sign(det(A)) = sign(det(i?))sign(det(P))sign(det(L))sign(det((5)). 

These signs of each these determinants is easy to compute. The norm of A~^ we 
compute with the power method. Essentially this method starts with a random nnit 





























































K-THEORY AND PSEUDOSPECTRA FOR TOPOLOGICAL INSULATORS 


39 




Figure 10.6. Phase transition in 3D. The left panel shows the disorder 
averaged new class All index. The right panel shows the disorder- 
average of gap localized at the origin. The number N of samples for 
an L-by-L-by-L lattice used in these averages was: L = 6, N = 26036; 
L = 12, N = 23357; L = 18, N = 3491; L = 24, N = 1547. 


vector V = Vo and then computes a few dozen iterations of 

1 

= Tl -n-Wn. 

l|Wn|| 


We modified this procedure a little. We found starting with vq having all entries 
equal worked better than a random vector in this setting. Then, following EH , we 
compute w„ using equation 


11.1 


and the Matlab \ operator that computes C 


- 1 , 


without inverting the matrix. Since 


{A*A)~^ = QU*L*PR-‘^P*LUQ*, 


and since {A*A) ^ and Q* [A*A) ^ Q, we can compute 

= V \ (f \ (P \ (R-^ \ {P‘ \{L\(U\ v„_i)))))). 


\-l 


Via 




Acknowledgments 

The author wish to thank Deborah Evans, Alexei Kitaev, Matthew Hastings, Joel 
Moore and Hermann Schulz-Baldes for illuminating discussions, mathematical and 
physical. 

This work was partially supported by a grant from the Simons Foundation (208723 
to Loring) and by financial support form the Erwin Schrodinger International Institute 
for Mathematical Physics. Most of the comupting was done on machines at the Center 
for Advance Research Computing at the University of New Mexico. 



















K-THEORY AND PSEUDOSPECTRA EOR TOPOLOGICAL INSULATORS 


40 


Reeerences 

[1] A. Altland and M.R. Zirnbauer. Nonstandard symmetry classes in mesoscopic normal¬ 
superconducting hybrid structures. Physical Review 55(2): 1142, 1997. 

[2] J. Bellissard, A. van Elst, and H. Schulz-Baldes. The noncommutative geometry of the quantum 
Hall effect. J. Math. Phys., 35(10):5373-5451, 1994. 

[3] David Berenstein and Eric Dzienkowski. Matrix embeddings on flat and the geometry of 
membranes. Physical Review D, 86(8):086001, 2012. 

[4] Jeffrey L. Boersema and Terry A. Loring. TC-theory for real (7*-algebras via unitary elements 
with symmetries. 

[5] Man Duen Choi. Almost commuting matrices need not be nearly commuting. Proc. Amer. 
Math. Soc., 102(3):529-533, 1988. 

[6] Timothy A. Davis and Iain S. Duff. An unsymmetric-pattern multifrontal method for sparse 
LU factorization. SIAM J. Matrix Anal. AppL, 18(1):140-158, 1997. 

[7] 1. S. Duff, N. I. M. Gould, J. K. Reid, J. A. Scott, and K. Turner. The factorization of sparse 
symmetric indefinite matrices. IMA J. Numer. Anal.^ 11 (2): 181-204, 1991. 

[8] Iain S. Duff. The design and use of a sparse direct solver for skew symmetric matrices. J. 
Comput. Appl. Math.., 226(1):50-54, 2009. 

[9] Spren Eilers and Terry A. Loring. Computing contingencies for stable relations. Internat. J. 
Math., 10(3):301-326, 1999. 

[10] Jean Bellissard Emil Prodan, Bryan Leung. The non-commutative n-th Chern number. 
arxiv:1305.2425. 

[11] A.M. Essin and JE Moore. Topological insulators beyond the Brillouin zone via Chern parity. 
Physical Review B, 76(16): 165307, 2007. 

[12] Andrew M. Essin and Victor Gurarie. Bulk-boundary correspondence of topological insulators 
from their respective green’s functions. Phys. Rev. B, 84:125132, Sep 2011. 

[13] P. Friis and M. Rordam. Almost commuting self-adjoint matrices-a short proof of Huaxin Lin’s 
theorem. Journal filr die reine und angewandte Mathematik (Crelles Journal), 1996(479):121- 
132, 1996. 

[14] I.C. Fulga, F. Hassler, and A.R. Akhmerov. Scattering theory of topological insulators and 
superconductors. Physical Review B, 85:165409, 2012. 

[15] I.C. Fulga, F. Hassler, A.R. Akhmerov, and C.W.J. Beenakker. Scattering formula for the 
topological quantum number of a disordered multimode wire. Physical Review B, 83(15): 155429, 
2011. 

[16] M. B. Hastings. Making almost commuting matrices commute. Communications in Mathemat¬ 
ical Physics, 291(2):321-345, 2009. 

[17] Matthew B. Hastings and Terry A. Loring. Almost commuting matrices, localized Wannier 
functions, and the quantum Hall effect. J. Math. Phys., 51(1):015214, 2010. 

[18] Matthew B. Hastings and Terry A. Loring. Topological insulators and (7*-algebras: Theory and 
numerical practice. Ann. Physics, 326(7):1699-1759, 2011. 

[19] Ilya Kachkovskiy and Yuri Safarov. On the distance to normal elements in (7*-algebras of real 
rank zero. arXiv preprint arXiv: 1403.2021, 2014. 

[20] CL Kane and EJ Mele. Z_2 topological order and the quantum spin hall effect. Physical review 
letters, 95(14):146802, 2005. 

[21] Vladimir V. Kisil. Mobius transformations and monogenic functional calculus. Electron. Res. 
Announc. Amer. Math. Soc., 2(l):26-33 (electronic), 1996. 

[22] Alexei Kitaev. Classification of topological insulators and superconductors. Lecture given at 
the IPMU Focus Week Condensed Matter Physics Meets High Energy Physics, University of 
Tokyo, 8-12 February 2010. 



K-THEORY AND PSEUDOSPECTRA EOR TOPOLOGICAL INSULATORS 


41 


[23] Alexei Kitaev. Periodic table for topological insulators and superconductors. In AIP Conference 
Proceedings^ volume 1134, page 22, 2009. 

[24] M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buhmann, L.W. Molenkamp, X.L. Qi, and 
S.C. Zhang. Quantum spin Hall insulator state in HgTe quantum wells. Science^ 318(5851):766, 
2007. 

[25] Markus Konig, Hartmut Buhmann, Laurens W. Molenkamp, Taylor Hughes, Chao-Xing Liu, 
Xiao-Liang Qi, and Shou-Cheng Zhang. The quantum spin hall effect: theory and experiment. 
Journal of the Physical Society of Japan^ 77(3), 2008. 

[26] Huaxin Lin. Almost commuting selfadjoint matrices and applications. In Operator algebras and 
their applications (Waterloo, ON, 1994/1995), volume 13 of Fields Inst. Commun., pages 193- 
233. Amer. Math. Soc., Providence, RI, 1997. 

[27] Terry A. Loring. Factorization of matrices of quaternions. Exposition. Math., 30(3):250-267, 

2012. 

[28] Terry A. Loring. Quantitative iC-theory related to spin Chern numbers. SIGMA Symmetry 
Integrability Geom. Methods AppL, lOiPaper 077, 25, 2014. 

[29] Terry A. Loring and Matthew B. Hastings. Disordered topological insulators via (7*-algebras. 
Europhys. Lett. EPL, 92:67004, 2010. 

[30] Terry A. Loring and Adam P. W. Sprensen. Almost commuting self-adjoint matrices — the real 
and self-dual cases, arxiv: 1012.3494. 

[31] Terry A. Loring and Adam P. W. Sprensen. Almost commuting orthogonal matrices. J. Math. 
Anal. AppL, 420(2):1051-1068, 2014. 

[32] Ian Mondragon-Shem, Taylor L. Hughes, Juntao Song, and Emil Prodan. Topological criticality 
in the chiral-symmetric aiii class at strong disorder. Phys. Rev. Lett., 113:046802, Jul 2014. 

[33] Yoshiko Ogata. Approximating macroscopic observables in quantum spin systems with com¬ 
muting matrices. Journal of Functional Analysis, 2013. 

[34] E. Prodan. Disordered topological insulators: a non-commutative geometry perspective. Journal 
of Physics A: Mathematical and Theoretical, 44:113001, 2011. 

[35] S. Ryu, A.P. Schnyder, A. Furusaki, and A.W.W. Ludwig. Topological insulators and super¬ 
conductors: tenfold way and dimensional hierarchy. New Journal of Physics, 12:065010, 2010. 

[36] Bjorn Sbierski and Piet W Brouwer. Z 2 phase diagram of three-dimensional disordered topo¬ 
logical insulators via a scattering matrix approach. Physical Review B, 89(15):155311, 2014. 

[37] Lloyd N. Trefethen and Mark Embree. Spectra and pseudo spectra. Princeton University Press, 
Princeton, NJ, 2005. The behavior of nonnormal matrices and operators. 

[38] M. Wimmer. Efficient numerical computation of the pfaffian for dense and banded skew- 
symmetric matrices. ACM Trans. Math. Software, 38(4), 2012. 

[39] Yan-Yang Zhang, Rui-Lin Chu, Fu-Chun Zhang, and Shun-Qing Shen. Localization and mobility 
gap in the topological anderson insulator. Phys. Rev. B, 85:035107, Jan 2012. 

Department of Mathematics and Statistics, University of New Mexico, Albuquerque, 

NM 87131, USA. 



